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Abstract 

For many biotechnological applications it would be useful to better under- 
stand the effects produced by electric fields on lipid membranes. This review 
discusses several aspects of the electrostatic properties of a planar lipid mem- 
brane with its surrounding electrolyte in a normal DC or AC electric field. 

In the planar geometry, the analysis of electrokinetic equations can be 
carried out quite far, allowing to characterize analytically the steady state 
and the dynamics of the charge accumulation in the Debye layers, which re- 
sults from the application of the electric field. For a conductive membrane 
in an applied DC electric field, we characterize the corrections to the elastic 
moduli, the appearance of a membrane undulation instability and the as- 
sociated flows which are built up near the membrane. For a membrane in 
an applied AC electric field, we analytically derive the impedance from the 
underlying electrokinetic equations. We discuss different relevant effects due 
to the membrane conductivity or due to the bulk diffusion coefficients of the 
ions. Of particular interest is the case where the membrane has selective 
conductivity for only one type of ion. These results, and future extensions 
thereof, should be useful for the interpretation of impedance spectroscopy 
data used to characterize e.g. ion channels embedded in planar bilayers. 

Keywords: lipid membrane, electric fields, electrokinetics, impedance 
spectroscopy 
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1. Introduction 

Bilayer membranes formed from phospholipid molecules are an essential 
component of the membranes of cells. The mechanical properties of equilib- 
rium membranes are characterized by two elastic moduli, the surface tension 
and the curvature modulus [l|, which typically depend on the electrostatic 
properties of the membranes and its surroundings Understanding how 
these properties are modified when the membrane is driven out of equilib- 
rium is a problem of considerable importance to the physics of living cells. 
A membrane can be driven out of equilibrium in many ways, for instance by 
ion concentration gradients or by electric fields. 

Quite generally one can distinguish between systems in which the electric 
field is applied externally and systems which are able to self-generate electric 
fields: 

1.1. Membranes in externally applied electric fields 

The external application of electric fields on lipid films is used to pro- 
duce artificial vesicles (electroformation), as well as to create holes in the 
membrane (electroporation) [sl] . Both processes are important for biotechno- 
logical applications and they are widely used experimentally. However, they 
are still not well understood theoretically. The research on electroformation 
is motivated by the hope to produce artificial lipid vesicles in a controlled 
and simple way, which will be key to many biotechnological applications. 
Cell electroporation is a popular technology and biomedical applications of 
in vivo cell electropermeabilization j3] are gaining momentum for drug and 
nucleic acids electrotransfer and for the destruction of tumor cells for cancer 
treatment jsj. 

In view of the importance of these applications, many research efforts have 
been devoted to study and understand deformations of giant unilamellar vesi- 
cles (GUVs) due to the application of electric fields. In the presence of an 
AC electric field, giant unilamellar vesicles show a rich panel of possible be- 
haviors and morphological transitions depending on experimental conditions 
- electric field frequency, conductivities of the medium and of the membrane, 
salt concentration, etc. 0,0]- A theoretical framework involving hydrody- 
namics and a continuum mechanics description of the membrane has been 
developed, which accounts quantitatively for the observed equilibrium and 
non-equilibrium shapes taken by the vesicles in the presence of an AC elec- 
tric field js], [oj . For a clear and self-contained presentation of this theoretical 
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framework, we recommend the chapter " Non-equihbrium dynamics of hpid 
membranes: deformation and stabihty in electric fields" by P. Vlahovska, in 
the same issue of this book. 

The application of external fields is also interesting as a means to move 



fluids via electro-osmosis [10|, IJJJ and to self-assemble colloidal particles, for 
various technological applications. Moreover, the ability to move fluids and 
nanoparticles at small scales is used in many biological systems. For instance, 
membrane-bound ion pumps and channels are able to transport water (for 
instance in aquaporin channels) and ions (in ionic pumps and channels) in a 
particularly selective and efficient way, which one would like to reproduce in 
artificial or biomimetic microfiuidic devices. 

1.2. Membranes in self- generated electric fields 

In some cases of biological relevance, membranes are able to self-generate 
an electric field, due to embedded ion channels or pumps. This can be 
achieved because the channels are able to transport ions from one side of the 
membrane to the other in a selective way, either down their concentration 
gradient in passive transport or against it in active transport, e.g. coupled 
to the hydrolysis of Adenosine triphosphate or activated by light. Probably 
the best known example is the opening and closing of ion channels in nerve 



cells allowing the transmission of an electric signal via action potentials [12 



For all these reasons, ion channels and pumps play an essential role in many 



biological functions of a cell [13 



In order to better understand how nerve cells operate in vivo, it would 
be helpful to construct an in vitro biomimetic equivalent which would have 
some key features of the in vivo system, such as the ability to generate an 
action potential, but without the complexity of a real nerve cell. Active 
membranes, which are giant unilamellar vesicles containing ion pumps such 



as bacteriorhodopsin [ij, ll5|, |16| are a promising system to achieve this goal. 



The main purpose of this review is to propose and analyze a simple model 
to foster the understanding of various effects resulting from electric fields 
acting on a planar lipid membrane. Although we are mostly interested in ap- 
plications to biological or biomimetic systems composed of lipid membranes, 
we would like to point out that the theoretical framework presented here is 
very general. It can be easily adapted to analyze the electrical properties of 
artificial membranes which can have very different properties from biological 
membranes (as far as e.g. ionic conductivities or the bending stiffness are 
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concerned) . 

This review is organized as follows: in section |2] we present the model 
for a planar lipid membrane and its surrounding fluid in an applied electric 
field. In this section, we will restrict ourselves to the case of a DC field. In 
particular we will focus on i) the electrostatic and electrokinetic steady-state 
corrections to the elastic moduli of the membrane due to the application of 
the electric field, see section 12. 3j ii) the flow fields which can be predicted 
from such an approach, at steady state and in the case that the membrane 
is ion-conductive, see section 12.41 In section 12.51 we will compare the model 
predictions to two relevant experiments. More details on this theoretical 
framework, as well as an extension to the nonlinear electrostatic regime us- 



17, 18 
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ing the Poisson-Boltzmann equation, can be found in Refs. 
Finally, in section |3] we present an analysis of the model in the presence of 
time-dependent AC electric fields. We provide derivations for the impedance 
of the system from the underlying electrokinetic equations, for situations 
where the membrane is either blocking or selectively conductive for ions. 



2. A quasi-planar membrane in a DC electric field 

The mechanical properties of membranes at equilibrium are characterized 
by two elastic moduli, the surface tension and the bending modulus. These 
moduli typically depend on electrostatic properties, and their modifications 
in the case of charged membranes or surfaces in an electrolyte have been 
examined theoretically in various situations: in the linearized Debye-Hiickel 
approximation as well as in the nonlinear Poisson-Boltzmann (PB) regime, 
for lipid monolayers and symmetric bilayers (2], [2l[ 22 , 23 1 . More recently, 
charged asymmetric bilayers with unequal Debye lengths on both sides of 
the membrane 2J] and an uncharged membrane in a DC field 25|] have been 
investigated. 

In all the works mentioned above, a free energy approach has been used. 
Note that while this method works well for equilibrium membranes, it is 
not applicable to situations in which the membrane fluctuations have a non- 
equilibrium origin, as in the case of active membranes containing ion channels 



26l 1271 . Il4l Il5| or in the case of a membrane in a time dependent electric 



field. In our recent work [17|, uM, uM l20|, we thus have studied this problem 



using an electrokinetic approach, which does not have the limitations of a 
free energy formulation. In this framework we allow for a finite conductivity 
of the membrane due to e.g. ion channels or pumps, and the ion transport is 
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described using a Poisson-Nernst-Planck approach [28|, |29|, |30[ . The electro- 
static corrections to the elastic moduli and the fluid flows in the electrolyte 
are then obtained by imposing the overall force balance at the membrane. 

Two additional points are worth emphasizing: first, our approach is able 
to correctly describe the capacitive effects of the membrane and of the De- 
bye layers while keeping the simplicity of the " zero-thickness approximation" 
on which most of the literature on lipid membranes is based. This is accom- 
plished by the use of an effective Robin-type boundary condition (BC) at the 
membrane. Second, as the method is based on a calculation of the general 
force balance at the membrane, additional non-equilibrium processes could 
be included into the model rather easily. For simplicity we investigate here 
only the effects of ionic currents flowing through the membrane, which in turn 
affect the fluid flow near the membrane. Other non-equilibrium effects that 
could be included as well are for instance including ion channel stochasticity 
or active pumping. 

2.1. Model equations - electrostatics 

Figure [1] shows a sketch of the planar geometry that is studied: we con- 
sider a steady current driven by a DC voltage drop V across two electrodes 
separated by a flxed distance L. The membrane is quasi-planar and located 
at 2; = 0. It is embedded in an electrolyte of monovalent ions with number 
densities and n~. It contains channels for both ion species but is itself 
neutral, i.e. does not carry flxed charges. The channels or pumps are as- 
sumed to be homogeneously distributed in the membrane and enter only in 
the effective conductance G, as introduced below. A point in the membrane 
is characterized in the Monge representation by the height function h{rj_), 
where is a two-dimensional in-plane vector. The base state of this prob- 
lem is a flat membrane. Hence the electric fleld, assumed to be perfectly 
aligned in ^-direction, is perpendicular to it. We assume a quasi-static ap- 



proach [25|, Il8| in which membrane fluctuations are much slower than the 
characteristic diffusion time r = for the ions to diffuse a Debye length. 
In the electrolyte, the electric potential obeys Poisson's equation 

1 2 
= — (en"*" — en") = — p. (1) 

Here e is the elementary charge, e is the dielectric constant of the electrolyte 
and we have introduced half of the charge density, 

72+ — n~ , . 

P = e ^ . (2) 
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Figure 1: Sketch of a quasi-planar membrane 
embedded in a symmetric electrolyte. The 
initially flat bilayer membrane is represented 
by the plane z = 0. The membrane fluctu- 
ations around this base state have not been 
represented. A voltage V is applied far from 
the membrane on electrodes separated by a 
distance L. The membrane carries ion chan- 
nels which give rise to a conductance G. 
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For the sake of simplicity, we assumed a symmetric 1 : 1 electrolyte, thus far 
away from the membrane = n~ = n*, and the total system is electrically 
neutral. The densities of the ion species obey the Poisson-Nernst-Planck 
(PNP) equations 

a^n^ + V • j± = , i^ = D Vn± T ^^-^^'p) ' (3) 

where are the particle current densities of the ions and fc^T is the thermal 
energy. We will assume here that both ion types have the same diffusion 
coefficient D. Note that we will discuss the effects of differing diffusion 
coefficients for an applied AC voltage in section 13.41 

Since we are primarily interested in the behavior close to the membrane, 
for the boundary conditions (BC) far away from the membrane we assume 

<P{z = ±L/2) = ±V/2, (4) 
p{z = ±L/2) = 0. (5) 

Eq. (jl]) states that the potential at the electrodes is held fixed externally. 
This BC is quite oversimplified for real electrodes, but captures the main 
effects of the electric field, see the discussion in Ref. 19|. We have also 
assumed that the distance between the electrodes is much larger than the 
Debye length, L ^ Xd = where 

S-^- (6) 
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Hence, as already mentioned above, the bulk electrolyte is quasi-neutral with 
negligible charge density (compared to the total salt concentration) and far 
from the membrane Eq. holds. 

The BC at the membrane is crucial to correctly account for capacitive 
effects. We use the Robin-type BC 



see 



Appendix A for a derivation) 



\m{n ■ V)0|2=/i+ = Am(n • V)0|^=/i- = 0(/i^) - (f){h ) 
where n is the unit vector normal to the membrane and 

Am d . 



(7) 



(8) 



Am is a length scale containing the membrane thickness d and the ratio of 
the dielectric constants e/e^ of the electrolyte and the membrane. Note 
that in Eq. (I7l), the membrane plays a similar role as the Stern layer in the 
description of Debye layers near a charged interface. This BC was rederived 
for electrodes sustaining Faradaic current 31, 33| or charging capacitively 
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18, 19 . There it was 



and was applied for membranes in Refs. [3( 
shown to properly account for the jump in the charge distribution which 
occurs near the membrane as a result of the dielectric mismatch between the 
membrane and the surrounding electrolyte. 

In addition to Eq. ([7]), we impose the continuity of the bulk current Jj^^o 
at the membrane. This BC involves the ohmic law 



4m 



=0 ; 



(9) 



where G denotes the membrane conductance per area and /i'' the electro- 
chemical potential. The electrostatic potential and the ion densities can now 
be obtained by solving Eq. ([1]) in the linear Debye-Hiickel approximation and 



one obtains 19 



i) the jump of the charge density at the membrane, p„ 

ii) the current through the membrane, jm, and 

iii) the electric field inside the membrane, EIP". 
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For simplicity, in the derivation of Eqs. f llOtfT^ we assumed equal ion con- 
ductivities {G^ = = G) and a symmetric electrolyte on both sides of the 
membrane (k^° = k,^^ = k). Note that the method presented in this section 
can be easily extended to cover more general cases. In additon, the nonlinear 
electrostatic problem (keeping the Poisson-Boltzmann equation) can be still 
solved analytically in the non-conductive case. The nonlinear generalizations 



of Eqs. can be found in Ref. |20 



2.2. Model equations - hydrodynamics and force balance at the membrane 

The hydrodynamics of the electrolyte is described by the incompressible 
Stokes equation, — Vp + ?7V^v + f = with V ■ v = 0, where v is the 
velocity field of the electrolyte, r] its viscosity, p the hydrostatic pressure and 
f = — 2pV0 the electric driving force. From the solution of the electrostatic 
and the hydrodynamic problem, one obtains the total stress tensor 



Tij 



-p6y + T] {diVj + djVi) + e (EiEj - ^SijE"^ , (13) 



which contains the pressure, the viscous stresses in the fluid and the Maxwell 
stresses. 

The lipid bilayer membrane, on the other hand, behaves as a two dimen- 
sional fluid which can store elastic energy in bending deformations. More 
precisely, its elastic properties can be described by the standard Helfrich free 
energy 

FH = \j ^^'rxpo [Vhf + Ko {y^h)\ (14) 

where Eq is the bare surface tension and Kq the bare bending modulus of 
the membrane. 

All forces present in the system, the electrostatic, viscous and elastic 
ones, have to fulfill the force balance equation. The latter states that the 
discontinuity of the normal-normal component of the stress tensor, as defined 
in Eq. ( !T3|) and evaluated at the membrane position, must equal the restoring 
force due to membrane's elasticity, hence 

dF 

- {rzz,i\z=h+ - rzz,i\z=h-) = "g^^ = {-^okl - Kokj_) h{k±) . (15) 

Here the index 1 in the stress tensor refers to the order of an expansion with 
respect to the membrane height field h{r±). Note that at zeroth order, the 
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membrane is flat and thus only electric forces and osmotic pressure balance. 
By expanding to linear order in the height field h{T±), and using 



h oc /ioe*''^"^+^^''^^* , (16) 

Eq. f|T5|) yields the growth rate s(k_L) of membrane fluctuations. Details of 
the derivation of s{k±) can be found in Refs. [l7|, ll8|, [l9|. We would like 
to emphasize that the force localized at the membrane surface is a priori 
unknown in this problem. Thus it must be determined self-consistently from 
the BCs for the velocity and the stress. 

2.3. Growth rate and renormalized elastic moduli 

The force balance Eq. (fT5|l determines the growth rate s{k±) entering the 
normal stress difference, 

rik^s{k^) = -^(So + AS)A;l-r«fci-i(iro + Air)A;i. (17) 

The electrostatic corrections to the surface tension, AS = AS^^ + AS^, and 
to the bending modulus, AK = AK,^ + AK^ can be decomposed into: 

i) an outside contribution due to the charges accumulated in the Debye layers 
and denoted with the index k; 

ii) an inside contribution due to the voltage drop at the membrane and 
denoted with an index m. They are given by 

for the contribution due to the Debye layers and by 

AS^ = -eUE^rd, = W - (19) 

for the contribution due to the field inside the membrane. 

Note that in Eq. ( ITTll . one also obtains a purely non-equilibrium correction 
Tk = It would correspond to a term proportional to in an "effective 

membrane free energy" incorporating the Maxwell stresses. At equilibrium 
such a term is forbidden by symmetry, but in a non-equilibrium situation, 
where the membrane sustains a current jm 0, it is allowed. For realistic 
parameters, however, this term is very small, see Ref. [18| for a detailed 
discussion. 
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Figure 2: The renormalized growth rate or 
dispersion relation, rs, as a function of the 
rescaled wave number k' = k^/n for three 
vohages: V = 0.7V (dashed hne), V = 0.75V 
(soUd hne), V = 0.8V (dash-dotted Hne). 
We have used the foUowing parameters: di- 
electric constants e = 80eo and = 2eo; 
membrane thickness d = 5nni leading to 
Am = j-d — 200nm; diffusion coefficient of 
ions D — lO^^m^s^-'^; viscosity rj — lO^^Pas; 
inverse Debye length k = 2 • lO^m"^; bare 
surface tension Sq = ImNrn"^; bare bending 
modulus K{) — lOksT. Here we assumed a 
non-conductive membrane, G = 0. 



The inside contribution to the membrane surface tension is always nega- 
tive, see Eq. f lT9|) . The same is typically true for the outside contribution, see 
Eq. ( lT8l) and note that Pm,jm > 0. Hence these contributions can overcome 
the bare surface tension Eq. If this is the case, an instability towards mem- 
brane undulations sets in. Such an instability had already been described 



for the high salt limit in Ref. [3J]. Note that the linearized theory devel- 
oped here describes only the early stages of the instability, but it is more 
general than previous works since it is not limited to the high salt limit and 
in addition accounts for hydrodynamic effects. The linear growth rate of the 
membrane fluctuations given by Eq. f|T7|) is shown in Fig. |2]in rescaled units. 
We scaled the wave vector by n, hence k' = k±/ k and the time by the typical 
time for ions to diffuse a Debye length, r = The control parameter of 
the instability is the external voltage V. Fig. [2] shows the growth rate for 
three different levels of the voltage: the dashed line is for V = 0.7V, which 
lies below the threshold of the instability, all wave numbers are damped and 
the membrane is stable. The solid and the dash-dotted line correspond to 
V = 0.75V and 0.8V. These values are above threshold and the growth rate 
is positive for a finite wave number window. 

For a more detailed discussion of the dependance of the corrections to the 
elastic moduli, the instability threshold and the characteristic wave number 
as a function of salt concentration and membrane conductivity, we refer the 



reader to Refs. 18, 19 
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Figure 3: Representation of the flows around 
the membrane beyond the instability thresh- 
old. The orientation of the electric field is 
towards negative values of z. Panel a) shows 
the flow generated by the membrane bending 
mode. Panel b) shows the ICEO flow. Fi- 
nally, panel c) shows the actual flow, which 
is the superposition of the former two and 
results in a strong flow near the membrane, 
oriented parallel to the surface. Both axes 
are scaled by the Debye length Param- 
eters are as in Fig. [2] except V — 3.165V, 
K = lO^m-i, G = lOS'm-2 and L = lO/im. 



2.J^-. Flow fields near a driven membrane 

We now summarize the main features of the fluid flows which arise near 
the membrane when it is driven by ionic currents 18(|. Fig. [3] was generated 



by selecting the fastest growing wave number and using the corresponding 
maximum growth rate. The shape of the membrane undulation is represented 
with the black solid curves. Fig. [3t) shows the flow field for a high membrane 
conductance and low salt, in the regime where the membrane is unstable 
due to the electrostatic correction to the surface tension and thus starts to 
undulate. The resulting flow is a superposition of two distinct flows: first, 
the typical flow associated to a membrane bending mode 35| as shown in 
Fig. E^). Second, the flow which results from the ion transport across the 
membrane. The latter flow has the typical counter-rotating vortices of an 



ICEO flow |36|, as shown in Fig. |3)d). Clearly, the superposition of these two 



flow contributions. Fig. results in a parallel flow close to the membrane, 
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in contrast to the usual bending mode flow given by Fig. 

For most realistic parameters - modest conductivities, not too low salt 
- the flow generated by membrane bending is usually dominating and hides 
the small ICEO contribution. To be able to observe the flows of Fig. [3l a high 
membrane conductance G and low salt are needed. Also, since for macro- 
scopic electrode distances L (e.g. of the order of millimeters), the voltage 
needed to induce the instability is very high, we have assumed a microscopic 
electrode distance of L = 10/im. While it might still be possible to observe 
flows for higher salt and macroscopic electrode separations, such situations 
can not be analyzed within the Debye-Hiickel approximation used here. 

Note that somewhat similar looking flow patterns have been recently ob- 
served experimentally in vesicles subject to AC electric fields in Ref. j6|. On 
closer inspection, however, it appears that these flows most probably have 
a different origin from the ICEO flows, since they are more likely to result 
from electrophoresis of charged lipids within the membrane. 

2.5. Applications to specific experiments 

Here we will briefly discuss how the framework presented above can be 
applied to recent experiments: the first experiment studied supported mem- 



branes subject to an electric field [37|, while the second one investigated 
active membranes IJ, llSl Il6 . 

S. Lecuyer et al. |37| recently performed neutron reflectivity measure- 
ments on a system consisting of two nearby membrane bilayers in an exter- 
nal AC electric field. One of the bilayers was close to the bottom electrode 
and used to protect the second one from interacting with the wall. The bare 
values of the elastic moduli were known from X-ray off-specular experiments 



for a similar system 38|, yielding Sq — 0.5mNm~^ and Kq ~ IbksT. The 
experiments were performed in an AC electric field at several frequencies. 
For the lowest frequency (lOHz) and for a voltage of V = 5V, the electro- 
static corrections to the surface tension and bending modulus were found to 
be AS ~ -1 ± 0.15mNm-i and AK ~ 185 ± IbknT. 

Assuming that the membrane is non-conductive, G = 0, and using an 
inverse Debye length of k = 2 • 10^m~^ (milli-Q water) and the experimental 
electrode distance of L = 1mm, our model yields AS ~ — 2-mNm~^ and 
AK ~ 190k bT. Thus the model successfully accounts for the order of mag- 
nitude of the electrostatic corrections observed in this experiment. Note, 
however, that the linearized Debye-Hiickel approach is not a good approxi- 
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mation in this applied voltages are rather high. For this reason, we 



recently extended our work to the Poisson-Boltzmann regime [20 . 

The second experiments we would like to discuss concerns active mem- 
branes, which are artificial lipid vesicles containing bacteriorhodopsin ionic 



pumps [IJ, |15|, |16| . These pumps are able to transfer protons unidirectionally 
across the membrane by undergoing light-activated conformational changes. 
The transport of protons across the membrane eventually builds up a trans- 



membrane potential. In Refs. [26|, |27|, ll5|, a hydrodynamic theory has been 



developed to describe the nonequilibrium fiuctuations of the membrane in- 
duced by the activity of the pumps. This work triggered substantial theo- 
retical interest in the problem, mainly focusing on the proper description for 
these non-equilibrium effects associated with protein conformational changes 



39, 40, 41, 42, 43 



In these models for active membranes, the electrostatic effects associated 
with the ion transport were not explicitly described. The framework pre- 
sented in this review provides a more detailed description of the ion trans- 
port, which could be useful to understand some aspects of active membrane 
experiments. From a contour analysis of giant activejvesicles, the fluctua- 
tion spectrum of the membrane was measured in Ref. [16], and a lowering of 
the membrane tension produced by the activity of the pumps was reported. 
Only the correction to the surface tension has been accurately measured in 
this experiment and many aspects of the transport of ions are still unknown. 
However, for simplicity we can assume that the passive state corresponds to 
a non- conductive membrane, G = 0, and the active state to a membrane 
with conductance G = lOSm"^. If we also assume a typical transmembrane 
potential of the order of 50mV, we can use the results for the corrections to 
the surface tension obtained above. Accounting for the rather high amount 
of salt using k ^ 5 • 10^m~^, we flnd a reasonable estimate for the observed 
tension lowering, AS ~ 3 • 10~^Nm~^. We also flnd that there is no mea- 
surable difference for the bending modulus between the active and passive 
state, in agreement with the experiments. The model further predicts a cur- 
rent density of jm — lAm"^ when the pumps are active, which corresponds 
to an overall current of IpA on a vesicle of size l/xm. 

This accord in orders of magnitude for the electrostatic corrections is 
quite promising. For a more detailed comparison between experiments and 
the presented model, it would be necessary to do experiments in varying 
conditions (ionic strength, conductance of the membrane, or orientation of 
the pumps in the membrane for instance). Combined measurements of the 
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membrane current and the transmembrane potential in the same experiment, 
using e.g. patch-clamp techniques, would also be desirable. 



3. Impedance of a planar membrane in an AC electric field 



Impedance spectroscopy [4^ is an effective tool to obtain a characteriza- 
tion of the electric properties of lipid bilayer membranes. The method has 
been used in particular for supported lipid bilayers, which are a promising ex- 
perimental system to characterize membrane proteins, channels or inclusions 
and more generally constitute the basis of highly sensitive detection tech- 



nologies, i.e. biosensors j45|. In the recent work 46 1, for instance, impedance 
spectroscopy has been used to characterize gramicidin D channels in pore 
suspending membranes. Nowadays, many biotechnology companies develop 
systems to measure the impedance of whole cells for e.g. screening or drug 
delivery. 

In many cases, the interpretation of the data obtained by impedance spec- 
troscopy is not that straightforward. Typically one uses equivalent circuits, 
which are sometimes controversial, since different models can be used for 
fitting the data. Moreover, there is often a lack of knowledge concerning the 
conditions of validity of these equivalent circuits to describe the diffuse charg- 
ing in electric Debye layers. To answer these questions, one possibility is to 
start with an electrokinetic description based on the Poisson-Nernst-Planck 



equations. With such an approach, the dynamics of diffuse charging [33 



and the current voltage relation in electrochemical thin films have been suc- 



cessfully analyzed [32|. This approach is also useful for relating impedance 
measurements to the properties of the diffuse layers near charge selective 
interfaces such as electrodes or ion-exchange membranes |47|. 

In the following, we extend the model studied in the previous sections to 
the case of an applied AC electric field. For simplicity the membrane will 
be assumed to be strictly planar and non- fluctuating. We use the Poisson- 
Nernst-Planck equations to evaluate the impedance of this system, which 
can be then compared to simple equivalent circuits. We will first present the 
generic time- dependent equations for the perturbation induced by the applied 
AC field. Then we proceed to calculate the impedance for the following cases: 
i) an ideally blocking membrane with equal diffusion coefficients for the two 
ion species, ii) the same system but with unequal ion diffusion coefficients 
and finally iii) an ideally non-blocking membrane which conducts selectively 
only one type of ion. 
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3.1. Time- dependent electric fields 

The Poisson-Nernst-Planck equations for an electrolyte have already been 
given in section 12.11 Taking the time derivative of the Poisson equation, 
Eq. ([1]), one obtains 

- eStVV = e{dtn+ - dtu-) = -e(V ■ j+ - V • j"), (20) 

where in the last equation, the conservation of ion densities, Eqs. ([3]), has 
been used. Through integration over space (assuming a one dimensional 
geometry), and using the definition of the electric field, E = — V0, it follows 



that 48, 49 



I = ed{E + eJ, (21) 

where the constant in of integration, I, is the total electric current density. 
The first term on the r.h.s. in Eq. fl2ip is the displacement current. The sec- 
ond term, J = j+ — = 2jp, is the particle current density. The displacement 
current was absent in the previous section because we assumed a stationary 
state, but for the time-dependent case it is crucial to obtain the response 
to an externally applied AC electric potential. We note that by virtue of 
the Poisson equation, Eq. ([T]), the total current density is divergence- free, 
V • I = 0, at all times. Further note that the experimentally measurable 
quantity is given by the total electric current. For this reason, it is the 
relevant quantity to calculate impedance as shown below. 

3.2. Equations for time-periodic perturbations of an equilibrium base state 

Let us assume an established equilibrium solution Cq{z), Cq{z) and 0o(^) 
for the electrolyte in the absence of the AC field, which could be caused by 
an additional DC field or a Nernst potential. For convenience we consider 
here the charge densities c^. Note that = en"^ and = 2eco/{ekBT). 
The equations for the electrostatic problem, see Eqs. ([H [3]) above, read 

ed'^^cj) = c- -c-^ , (22) 
dtc^ = -d.J^, (23) 

= • (24) 

Linearization around the base state like 

C+ = 4 + T]C+ , c" = Co + VC~ , = 00 + , 
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where t] is a small book-keeping parameter, leads at order 0{ri^) to 

e(/)Q (z) 

= coe"^^^, with 0o solution of ed^cpo = cq (e<^° - e-"^") . 

This restates that the equilibrium solution has to fulfill the classical PB 
equation. At order 0{r]^) in the perturbations, we get 

e5> = C--C+, (25) 



d^C^ = D^d,[d,C^Tctj^d,^TC^j^d,<P,j . (26) 

As already discussed in the general case above, taking the time-derivative of 
Eq. (I25ll . insertion of the linearized PNP equations ( l26l) and integration in z 
yields 

+ (^d.C+ + 4-^0.^ + C+-^a,(/.o) = lit) . 

The integration constant I{t) is the total electric current density. As we 
are interested in the response to an AC external voltage, V{t) = Vqc*'^*, 
introducing I{t) = /qC*^* and $ cx e*'^*, we arrive at 

+D+d,C+ - D-d,C- + + D-C-) y^9,0o = h ■ (27) 

The first term on the l.h.s. is the displacement current. The remaining terms 
are currents due to concentration gradients and a current induced by the equi- 
librium potential at the membrane. All these contributions taken together 
yield the total current Jq in response to the external AC field. 

We are left with the problem to solve Eqs. f l26|) and fl27j) with the external 
voltage entering via the boundary conditions, just like in section |2j 

3.3. Impedance for an ideally blocking non- conductive membrane 

The equations derived in the last section are general as they describe the 
first order perturbation in an electrolyte induced by an AC voltage externally 
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imposed at some boundaries. Let us now apply them to the planar membrane 
geometry as sketched in Fig. [H The membrane is assumed to be flat and 
located at 2; = 0. The AC voltage will be externally applied at the electrodes 
aX z = ±L/2. For simplicity, we assume that there is no additional DC 
electric field or Nernst potential, i.e. that the equilibrium solution is given 
by the homogeneous solution 0o = 0, Cg = Cq. 

First we will treat the simplest case of an ideally blocking, non- conductive 
membrane, j^(0) = 0. We also assume equal diffusion coefficients for the 
positive and negative ions, = = D. Then the above equations (126| 
[27|l for the perturbations reduce to 

tuC^ = Dd, (^d,C^ T ^o^d,^^ , (28) 

icoe + 2^Co-^^ + D{d,C+ - ^.C") = /q . (29) 

Due to the symmetry of our system, one has 

<l>{z,t) = -<^{-z,t) , p{z,t) = -p{~z,t) , c{z,t) = c{-z,t) . (30) 

Hence it is enough to solve the problem in z E [—L/2,0]. The BCs in the 
chosen geometry read 

C+(-L/2) = 0, (31) 

C-(-L/2) = 0, (32) 

$(-L/2) = -K,/2, (33) 

d^C^iO) + coj^dMO) = = j+(0)/D , (34) 

d,C-iO) - coj^dMO) = = riO)/D , (35) 

X^mdMO) = $(0+) - H0-) . (36) 

Eqs. fl3m33p fix the densities and the potential at the electrodes, as has 
already been discussed in section 12.11 The next two equations flMj |35|) state 
that the membrane is non-conductive for both ion types. Finally the last 
equation f l5B]) is again the Robin-type BC describing the capacitive behavior 
of the membrane with the effective length scale = —d. We will use the 
first five BCs to fix the five integration constants of Eqs. ( !28l [29]) . Then 
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imposing the last condition will yield the current-voltage relation and finally 
the impedance. 

Extracting an equation for Cs = + by adding Eqs. f l28|) yields 
iujCs = DdlCs. From the BCs d.Csifi) = = C,(-l/2) it follows = 0, 

i.e. the total density of particles (positively and negatively charged) remains 
homogeneous. Introducing p = — and substracting Eqs. f l25]) yields 

iujp = Ddlp + DtK^dl^ , (37) 
{iue + DeK^) d,^ + Dd,p = Jq , (38) 

where we have used = et^ 12. Eq. ( l38i) can be integrated, yielding 

DtK^ + luje 

The BCs p{-L/2) = 0, <l>(-L/2) = -Vo/2 fix the constant of integration to 

Vo ^ IoL/2 



2 Dsk"^ + ioje 
Insertion of the obtained potential into Eq. fl37|) for p yields 

^ P = '9fp. (39) 

Using again the obtained potential transforms the EC dzp(fi) + eK^dz^{Gi) = 
into the simpler form dzp{0) = i^^. Together with p(—L/2) = 0, the 
solution of Eq. (l39|l can be given as 

with the (complex) inverse length scale 



/3 = ^/t2 + ^cu/D . (41) 



The remaining BC, Eq. (1361) . is a jump condition at the membrane. What 
we have calculated above are the solutions on z G [— L/2,0]. Using 

the symmetry of our problem, Eq. fl30l) . one directly obtains on 
z G [0,L/2]. Imposing Eq. ([36]), A^9^$(0) = <l>>°(0) - <l><°(0), then yields 



DdzPiO) _ .0 , -Dp>'{0) 



\ — -^^rv-y _ >u _|_ 
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.0 , -^P<°(0) 



Solving for the external voltage Vq - note that it enters in the integration 
constant Ci of the electric potential - one gets 

,r _ hL , Dp{^) , , Jo - Dd-,p{Q) 

DtK'^ + lUJt DtK'^ + loot DtK'^ + lUJt 

This is the current-voltage relation. The impedance is defined as Z{u:) = 
V{u)/AI{u) = Vo/{AIo), with A the membrane area normal to the 2;-direction. 
Using the expression for the density, Eq. fHU]) . one arrives at the following 
expression for the impedance of a non-conductive membrane 

L/A . (L/A)^ tanh[/j(L/2)] {XJA)(l-^f^) 



Let us discuss the obtained result. The first term is the contribution of 
the electrolyte. This can be seen by rewriting it as 

1 



(43) 



and identifying the capacitance of the bulk, Cb = ^A/L, which is in parallel 
with the resistance of the bulk Rb = 75^^ = 2DcoeA ■ ^ similar interpreta- 
tion holds for the term j^^^^^^^ in Eq. fj42l) . which can be written as 

Zs = „-i !_ . ^ . (44) 

This is again a RC-circuit with the capacitance Cs = eA/X^ = ^mA/d of 
the membrane and a resistance Rs = j^^- It arises from the Robin-BC 
which involves the effective length scale defined in Eq. ([H]). One can thus 
recast Eq. f H2|) into the form 



7 +7 ' 

+ ^TITf 



uj Rn + iujCb 



tanh[/3(L/2)] 

I3L/2 L 



(45) 



The last term in this equation, let us call it Zc, is due to charging of the 
double layer and the membrane. This can be best seen in the limit Ad/L = 
l/(/tL) ^ 1, i.e. when the Debye length is small compared to the system 
size. Then ^^'^^Wj^J'^)^ ^ 2/{kL) and in the prefactor, the resistance R'^^ 
dominates over the capacitance. One gets 

Zc ^ — ^ , (46) 

lUjCeff 
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Figure 4: Effective circuit for the ideally blocking non-conductive membrane, Eq. (j47|): 
Two RC-circuits, one for the bulk and one for the membrane are in series with the effective 
charging capacitance of the membrane. 



with the effective capacitance C^ff = e 2\J\^x ■ Note that the thickness of the 
corresponding planar capacitor is the sum of the two Debye layers thicknesses 
(2 A/)) and the effective length describing the capacitive effects of the 
membrane. 

As shown in Fig. HI for the blocking non-conductive membrane one effec- 
tively has an association in series of the RC-circuit of the bulk, the RC-circuit 
of the membrane and the effective capacitance of the charging membrane 

Z = Zb + Zs + Zc= _ 1 ,\ ^ + „_i ^ . „ + , (47) 

as long as Xd/L <C 1 holds. As Am — 200nm, the impedance contribution 
Zs is usually small compared to the bulk resistance and can be neglected 
for L ^ Xm- However, the contribution described by Am to the charging 
impedance Zc can be of similar order as the one from the Debye layers and 
might even dominate the charging. 

The best way to visually present the impedance is by a so-called Nyquist 
plot |4^]. There one traces the negative imaginary part, — Im[Z(a;)], of the 
impedance as a function of its real part, Re[Z(a;)], for varying frequency u. 
Nyquist plots for the full impedance, Eq. ( H5|) . and for the limit Xd/L <C 1, 
Eq. ( H71) , are shown in Fig. [5l Panel a) shows the case of a macroscopic system 
size, L = 1mm. One clearly notices the RC- semi- circle terminating for high 
frequencies at the origin. For the given parameters one enters this semi-circle 
at a; ^ 50Hz; the maximum is achieved for urc = Dk^ = IkHz. The low 
frequency branch is dominated by the membrane charging capacitively at 
R ~ Rb + Rs, thus for low frequencies one has a divergence like {iujCeff)~^ ■ 
As Xd/L ~ 10~^, the effective circuit and the full calculation agree well. 
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Figure 5: Panel a) shows a Nyquist plot for a macroscopic system size, L = 1mm. At 
low frequency the behavior is governed by the charging of the membrane. The semi-circle 
is governed by the bulk-RC-circuit. As \d/L ~ 10~^, the effective circuit and the full 
calculation agree well. Panel b) shows a Nyquist plot for a microscopic system, L — 10/zm. 
As L decreases, the bulk becomes less important and the RC semi-circle less pronounced. 
The full calculation (solid curve) yields a lower resistance for the charging process at low 
frequencies than the effective circuit (dashed curve) . Parameters as in Fig. [2] except for 
K = 10~^m~^ (pure water); membrane area A — IfiTc? . 



Fig. [5b) shows the case of a microscopic system size, L = lO/im. Here the 
bulk RC-signal is much less pronounced and charging dominates entirely. The 
full calculation (solid curve) yields a lower resistance for the charging process 
at low frequencies than the effective circuit obtained by the small-Debye layer 
approximation (dashed curve). 

3.4- Non- conductive membrane: effect of unequal diffusion coefficients 

In this section we investigate the effect of differing diffusion coefficients 
for the two ion species, ^ , on the impedance of a blocking non- 
conductive membrane. Except for this assumption, the calculation is analo- 
gous to the one of the previous section. Equations ( l26i) . ( 1271) for the pertur- 
bations now read 

(48) 

D-d,C- = Jo . (49) 

The BCs are still given by Eqs. ( l3Tll36ll . Since the equations for the charge 
densities do not decouple as before, it is useful to introduce C = + C~ 



en 
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and p = C+ — C again, yielding 

iooC = dl {DC + 5p + 25c<l>] , 
%up = dl {bC + Dp + 2Dc^] , 
{iue + 2Dc) 9,$ + 6d,C + Dd,p = Jq . 

Here we introduced the average and the difference of the two diffusion coef- 
ficients 

D = {D+ + D-)/2, 6 = {D+ - D-)/2. (50) 
Integration of the equation for the potential $ yields 

^ ^ /o.-(^CM+gpM) ^ /«i/2 



Den'^ + iwe 2 De/t^ + iue 

Insertion into the equations for C and p yields a matrix equation 

f ^UJ-[D-6W/N]^l -[6^UJe/N]^! \ f C\ _ 

\ - [6tue/N] dl tu - [Dtue/N] dl J ' \ p J ' ^ ' 

where we have introduced = DeK^ + iue. 

Assuming solutions of the form C,p oc e^^ , Eq. f l5T]) yields four solutions 
for the decay length /3. In the case of equal diffusion coefficients studied 
previously, (5 = and the equations are decoupled. In that case D = D and 
one simply gets (3f = ^ associated to the relaxation of the total particle 
density C and 1^1 = ^ = '^''^"^ associated to the relaxation of p, see 
Eq. ( 15^ . In the case of unequal diffusion coefficients, the equations are 
coupled and the general solutions are 



tuD + {D^ - 6^) K^I2 T y^(/€72)2(D2- 52)2 _ ^2^2 

= (D2 - 52) • (52) 

Here the minus sign applies to f3i and the plus sign to (32. Consequently, 
Eq. (ISTll is solved by the ansatz 



C 
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with the eigenvectors given by 



[6icoe/N] f3f 



[D - d^enyN] /3| 



The effective BCs read: d,C{0) = and d,p{0) = at z = 0; C{-L/2) = 
and p{-L/2) at z = -L/2. The last two BCs yield E^Bi + E2B2 = 
and Bi + B2 = 0. As Ei ^ E2 this implies Bi = = B2, i.e. the cosh- 
contributions in the solution vanish. After some algebra one obtains 



P 
C 



./o/t' E,E2 



sinh [f3i {z + L/2)] _ sinh {^2 {z + L/2)] 
u - El [ Ei/3i cosh(/3iL/2) cosh(/32L/2) 



.IqK^ E1E2 



u) E2 — El 



sinh {z + L/2)] sinh [/32 {z + L/2)] 
/3i cosh(/3iL/2) /32 cosh(/32L/2) 



Using the Robin-type EC, Eq. (136|1 . and once again the symmetry of the 
problem one gets 



A. 



{5d,C{0) + Dd^pjO)) 
2Dc + iuje 



-2cf - 2 



- (5C<0(0) + 5p<0(0)) 



2Dc + iue 



Solving for Vq, insertion of the obtained solutions for C and p and applying 
Z = Vq/^IqA) one obtains the impedance 



L/A i>^rn/A) (1 - DZ^ 



2/tV^ 



— t 



+ i 



Den? + icoe 
E1E2 



{Dk^ + iuj)uje E2 - El 

2k l;iL;2 

(L)/€2 + E2 - El 





^ tanh[/3iL/2] 






tanh [(32L/2] 


/32 



(53) 



The first two contributions are already familiar to us, they stem from the 
bulk and the Stern-like description of the membrane. Note that D enters 
instead of D. 

Let us discuss the newly arising terms. As an expansion in A^/L ^ 1 is 
a bit tedious, let us consider only the simpler limit k — > 00. Eq. ( 152|) for f3f 
has a minus sign in front of the square root, the two K-terms cancel and 



iujD 
(L)2 - 52) 



^ (5i = ±Jiuj/D^ 



(54) 
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Figure 6: Effective circuit for the ideally blocking non-conductive membrane with differing 
diffusion coefficients, Eq. ((55|) : Two RC-circuits, one for the bulk and one for the membrane 
are in series with the effective charging capacitance and a Warburg resistance. 



with Deff = (-D^ — 5^) /D. For /3| one has the plus sign in front of the 
square root, the terms in dominate and one simply gets (32 = ±k. For 

the eigenvectors to leading order one has Ei = - — — , E2 = and 



E2-E1 y ^ El J k2(s2-<52) ' E2-Ei V y 

Consequently, the last term in Eq. ( I53l) exactly reduces to the Debye-layer 
part of the charging contribution. Finally one obtains at leading order in Ad 



with 



■ = Zb + Zs + Zc + Zw (55) 
2X1 tanh [PiL/2] 



DeA/L [D^ - 52) p^L 



(56) 



The first two terms are the RC-contributions of the bulk and the membrane 
(note that D = {D^ + D^)/2 enters instead of D). The charging capacitance 
Zc of the membrane is unchanged. The last term is the so-called Warburg 
impedance, with /3i = ^/ioo/ D^f / and -De// = (-D^ — 5^) / D. Note that this 
contribution is only present for unequal diffusion coefficients 5 = {D~^ — 
D^)/2 7^ 0. It is proportional to A|, at leading order 0. 

The effective circuit corresponding to the obtained impedance is shown in 
Fig. [6l The contribution Zyi/ has been first described by Warburg 5^, 51] for 



"'^For simplicity, we used the limit k ^ 00 to derive this term. Taking this limit strictly, 
the contribution would vanish ~ as then both charge types diffuse infinitely rapidly across 
the zero-thickness Debye-layer. In real systems, however, k remains always finite and thus 
one should include the leading order contribution, Zw, in the impedance. 
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electrochemical systems; in a nutshell, it arises from damped concentration 
oscillations close to an interface, here the membrane. We note however, that 
with typical differences in diffusion coefficients / D~ = 0.1. .10, a Nyquist 
plot of Eq. fl55|) is indistinguishable from Fig. |5] obtained for equal diffusion 
coefficients. This is due to the fact that in the geometry under investigation, 
the charging of the membrane is highly dominating the low-frequency be- 
havior as it is proportional to u'^. Nevertheless, experiments often display 



a Warburg-like impedance at low frequencies, see e.g. Ref. [46|. In the next 
section we will investigate the case of a slightly conductive ion-selective mem- 
brane and will find that in this case one indeed obtains a Warburg impedance. 
We thus postpone a discussion of to the next section. 

3.5. Impedance for an ideally non-blocking conductive membrane 

For many applications it is interesting to account for a small but non- 
zero membrane conductivity. This is important for instance in the context 
of the characterization of ion channel proteins or pumps embedded in a lipid 
membrane using impedance spectroscopy. In contrast to section 12.11 where 
we discussed the effects of a DC voltage on a conductive membrane that lets 
pass both types of charged ions (G"*" = G~ = G), here we will treat the case of 
a selective membrane, which lets pass only the positive ions. Thus, we assume 
a linearized relation j+ = G^Afx^ where G~^ is the effective conductance per 
unit area. The negative ions are not allowed to pass the membrane, hence 
j~ = or effectively G~ = 0. This situation is relevant for biomembranes, 
where ion channels allow the passage of positively charged ions like Na"*" or 
K"*", but not of negatively charged ions like Cl~ which are typically larger. 

To simplify the analysis, we will not describe the structure of the Debye 
layers as explicitly as in the previous sections. Instead we rely on two known 
approximations used in the study of electrochemical systems: 

i) the bulk is to a good approximation locally electroneutral. More pre- 
cisely, deviations from electroneutrality occur only in the third order in an 
expansion of Xd/L, which is very small for usual system sizes. This result 
can be obtained using a matched asymptotics expansion [33| . Consequently, 
we will assume for all z, p{z) = 0, or C{z) = C'^{z) = C~{z) for the pertur- 
bation of the charge densities. 

ii) Although we do not treat the Debye layers explicitly, we still im- 
pose effective boundary conditions for the electrochemical potential at the 
membrane. Thus, we implicitly assume that the electrochemical potential is 
continuous across the Debye layers. 
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We keep the geometry as before, i.e. a flat membrane located at z = 
with given AC voltage Vq at the electrodes located at 2; = ±L/2. We 
again assume that there is no additional DC field or Nernst potential, and 
equal diffusion coefficient^ for positive and negative ions. Using the above- 
discussed approximations, we obtain 

iuC = DdlC , (57) 
[iue + Dek^) d,^ = Iq . (58) 

Eq. f l58|) is again easily integrated for z G [— L/2,0] and together with the 
BC <l>(-L/2) = -1^0/2 one gets 

tooe + Dek'^ 

In addition we need three more BCs, namely 

C{-L/2) = 0, (59) 
D(dMO)-cdMO)) = J- = 0, (60) 

Did.CiO) + cdMO)) = j+ = — f^[C]o + e[$]oV (61) 

where [C]o = C{0^) — C(0~) and analogously for [$]o. The second condition 
is the no-flux condition for the anions. The third condition states that the 
bulk current of cations equals the current through the membrane, and is 
assumed to follow Ohm's law. From Eqs. (157 1 159 1 160|) . we obtain the following 
frequency dependent ion density distribution 

C(z) = — ; — — r- — — — -— sinh (a(z + L/2)) , 

^ ^ 2a{iuje + Den^) cosh {aL/2) ^ ^ ' 



where a = ^JiujjD is of Warburg-type, cf. Eq. ( 15^ . Note that here the 
Warburg impedance arises from breaking the cation/anion symmetry, due to 
differences in membrane conductivities rather than due to differences in their 
diffusion coefficients as in the previous section. Also note that although the 



^ Note that in case of unequal diffusion coefficients, one gets a contribution like dzC 
in Eq. (|58|). The subsequent calculations can still be performed in a completely analogous 
way. 
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Figure 7: Effective circuit for the ideally blocking and selectively conductive membrane, 
Eq. (|62l) : Two RC-circuits, one for the bulk and one for the membrane are in series with 
a Warburg resistance, caused by the ion selectivity of the membrane. 



membrane is non-conductive for the anions, this is a collective effect in which 
both types of moving charges participate. 

Finally, we use the BC for the cationic current, Eq. fl6T|) . to solve for the 
voltage Vo and obtain the impedance via Z = Vo/{IoA) as before 

, DKy{G+A) , fcBT/.V(ecoA)tanh(aL/2) 

^ - + ^D^^^T^ + ^^jt^d;;;^ a — • ^^^^ 

Here we already have identified the bulk circuit, it is present as in the previous 
cases. The second term is the membrane contribution. It can be written as 

Zm = p-1 \ ^ , (63) 

with the membrane's resistance, Rm = ^/{G~^A), and capacitance. Cm = 
The third term is the Warburg impedance, reading 



'W 



2\j^ tanh {^\Jiuj / DL/2 
DtA/L ^iuj/DL 



(64) 



for small uj. Note that it is of the same form as Eq. ( !56l) obtained for unequal 
diffusion coefficients, except for that in the latter appears an additional factor 
containing the diffusion coefficients. 

We can conclude that as a result of the ionic membrane selectivity, a 
Warburg impedance arises. Fig. [7] shows the effective circuit. A Nyquist plot 
is given in Fig. [HI One can identify the typical shape of a Warburg impdeance 



27 



l.Oxlo'^ 




Figure 8: Nyquist plot for a selectively conducting membrane. At high frequencies one has 
an RC-semi-circle, which is either dominated by the bulk or by the membrane, depending 
on the membrane conductance and the dimensions of the system. The low frequency 
behavior is governed by the Warburg impedance. Parameters as for Fig. [5] except for: 
L = 1mm; k = 10~^m~^ (pure water); A = Ijiw?. 

for low frequencies: namely, for decreasing frequencies, on leaving the RC- 
signal of the bulk — Im[Z(a;)]/Re[Z(a;)] acquires a slope of 45°. Finally, due 
to the finite system size Im[Z(ci;)] vanishes for — )■ 0. 

As already stated above, the calculation in this section [375] is oversimpli- 
fied. By assuming that the electrochemical potential is continuous across the 
Debye layers, there is no explicit contribution from the charging of the De- 
bye layers to the impedance. Hence A^, which is important for the charging, 
does not enter - indeed we did not even use the Robin-type condition. As the 
membrane is conductive, at least for the cations, charging of the Debye layers 
is of minor importance for the overall impedance. With a proper treatment 
of the charging of the Debye layers, using a matched asymptotics calculation, 
the Robin-type condition will reoccur to match the two solutions and will 
reintroduce the length scale into the problem. 

4. Conclusion 

The study and theoretical description of the effects induced by electric 
fields on lipid membranes in an electrolyte is a vast, challenging and far 
from fully explored problem, which is of relevance for many applications in 
biotechnology. 



28 



In this review, we have presented a theoretical framework to understand 
some of these effects in the simple case of a planar geometry. We have seen the 
importance of capacitive effects, occurring as a result of charge accumulation 
in the vicinity of the membrane, leading to renormalized elastic moduli and 
to membrane instabilities. We also have analyzed the flow fields which can 
be induced by currents due to small membrane conductivities. We discussed 
these effects only for time-independent (DC) electric fields, i.e. in the steady- 
state regime. 

Clearly, time-dependent electric fields lead to capacitive charging of the 
membrane and to time dependent membrane dynamics. The capacitive 
charging can be described in two ways: the first approach is based on the 
leaky dielectric model developed by G. I. Taylor [52j]. This approach is ex- 
plained and illustrated in the contribution of P. M. Vlahovska in the same 
issue of this book. One advantage of such an approach is that it captures the 
main physical effects associated with capacitive charging without the com- 
plexity of models which deal explicitly with the ion concentration fields. For 
this reason, it is useful to describe for instance the complex shape changes 
occurring in closed lipid vesicles 

The second approach, which we used in this work, is based on the elec- 
trokinetic Poisson-Nernst-Planck equations. This more refined level of de- 
scription includes ion concentration fields, and therefore it is useful to de- 
scribe specific effects associated for instance with the ion transport in ion 
channels or for effects occurring in low salt conditions. It is also needed 
to describe more precisely the capacitive charging, which as we have shown 
here includes several contributions coming from the bulk, the membrane 
impedance and the Debye layers themselves. In this review, we have tried 
to illustrate the strength of this approach for quantifying the impedance of 
a membrane-electrolyte system. In particular, we have shown how effective 
circuits used to interpret experimental data can be directly derived by this 
method. The membrane selectivity with respect to ion species is crucial 
to understand the conduction properties of membranes with embedded ion 
channels. We hope that our work will motivate further experimental and 
theoretical investigations in this field. 
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Appendix A. Robin-type boundary condition 

In brief, this boundary condition can be motivated for a fiat membrane as 
follows: since the membrane is assumed to bear no fixed charges, the normal 
components of the electric displacement are continuous at the two interfaces 
between the membrane and the electrolyte, 

ed,(t){z = ±d/2) = emd,(t)m{z = ±d/2) , (A.l) 

where (pm is the electric potential inside the membrane. Since the electric 
field Em = —dz(f)m is constant (to leading order) inside the membrane, the 
integral of the inside field can be written in the following way 

I'd/2 

Emd = / E„4z = - [Md/2) - (pm{~d/2)] = - [^{d/2) - 0(-rf/2)] , 

J~d/2 

where in the last step we used the continuity of the potential at the membrane 
surface. Together with Eq. (lA.ip this yields 

Xmd.^iz = ±d/2) = 0(rf/2) - (pi-d/2) . (A.2) 

If we take the limit — > while keeping = —d constant, one obtains 
Eq. d?]) in the particular case of /i = and n = z. The same derivation holds 
for the case of a slightly perturbed membrane surface h{r^), where is the 
in-plane vector. 
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